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Abstract. We study bifurcations of periodic orbits in three parameter general un¬ 
foldings of certain types quadratic homoclinic tangencies to saddle fixed points. We apply 
the rescaling technique to first return (Poincare) maps and show that the rescaled maps 
can be brought to a map asymptotically close to the 3D Henon map x = y,y = z,z = 
Ml + M 2 y + Bx — z"^ which, as known [1], exhibits wild hyperbolic Lorenz-like attractors 
in some open domains of the parameters. Based on this, we prove the existence of inhnite 
cascades of Lorenz-like attractor^. 
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1 Introduction 

In [1] it was discovered that the three-dimensional Henon map 

X = y, y = z, z = Ml + M 2 y + Bx — z^, (1.1) 

where [Mi, M 2 , B) are parameters {B is the Jacobian of map), can possess strange at¬ 
tractors that seem very similar to the Lorenz attractors, see hg. [1] Later it was shown 
that such discrete Lorenz attractors can arise as result of simple, universal and natural bi¬ 
furcation scenarios realizing in one-parameter families of three-dimensional maps IHEI]. 
This means, in fact, that the discrete Lorenz attractors can be met widely in applications. 
For instance, in [521153] such attractors were found in nonholonomic models of rattleback 
(called also as a Geltic stone). See also [411153] where various types of strange homoclinic 
attractors, including discrete Lorenz ones, were investigated. 

The similarity between the discrete and classical Lorenz attractors appears to be not 
accidental and it can be explained by various reasons. Thus, it is well known that the 
classical Lorenz attractor can be born as a result of local bifurcations of an equilibrium 
state with three zero eigenvalues when a flow possesses a (Lorenzian) symmetry [45] . 

^This work was supported by grant 14-41-00044 of the RSF and grant of RFBR 13-01-00589. 
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Figure 1: Plots of attractors of map (jl.ll) observed numerically in [I] for Mi = 0,5 = 0.7 and 
M 2 = 0.85 (left) or M 2 = 0.815 (right). In the left panel, the projection on the (a;, ?/)-plane is also 
displayed. In the right panel, a ’’figure-eight” saddle closed invariant curve inside the lacuna is shown. 
Note the similarity to the Lorenz attractors of the Shimizu-Morioka system, see HH [551156]. 


Analogously for maps, discrete Lorenz attractors can arise under bifurcations of fixed 
points with multipliers (—1,—1,-M), in this case the required local symmetry exists au¬ 
tomatically due to negative multipliers. As it was shown in mm, the second iteration of 
the map near this point can be embedded into a flow up to asymptotically small periodic 
non-autonomous terms. The corresponding flow normal form of such bifurcations coin¬ 
cides with the well-known Shimizu-Morioka model, which, in turn, exhibits the Lorenz 
attractor for certain parameter values [311EO]- Thus, we can consider the attractor in the 
map as the one of the Poincare map (period map) of a periodically perturbed system with 
the Lorenz attractor. On the other hand, as it was shown in paper |3^ by Turaev and 
Shilnikov, such discrete attractor is genuine in the sense that every its orbit has positive 
maximal Lyapunov exponent^. 


^Moreover, it is a wild pseudohyperbolic attractor (42], since it allows homoclinic tangencies (i.e. 
contains Newhouse wild hyperbolic sets [49]) and has an adsorbing domain inside which the differential 
DT, for the map itself and all close maps can be decomposed into a direct sum of transverse invariant 
subspaces W®'* and where is strongly contacting and expands exponentially volume. 
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In the present paper we study bifurcations of three-dimensional diffeomorphisms with 
homoclinic tangencies, leading to the birth of descrete Lorenz attractors. Problems of 
this kind were previously analyzed in [l3l H6l Wf[ HSl 0]. 

In 071 08] the birth of discrete Lorenz attractors from nontransversal heteroclinic 
cycles of three-dimensional diffeomorphisms was studied. Such a cycle contains two hxed 
points Oi and O 2 of type (2,1), i.e. with dimfP^(Oj) = 2,dimIP“(Oj) = 1, and one pair 
of stable and unstable manifolds intersect transversely and another pair has a quadratic 
tangency. It was assumed that at least one of the points Oi and O 2 is a saddle-focus, 
see Fig. Moreover, in all cases the additional condition that the Jacobians of the map in 
points Oi and O 2 are greater and less than one respectively was imposed (the so-called 
case of contracting-expanding maps). The birth of discrete Lorenz attractors was proved 
for three-parameter general unfoldings. 

Remark. Naturally, three parameters are needed to allow generically the existence 
of triply degenerate hxed points in the corresponding hrst return maps. In such families 
the hrst return map can be rescaled to the form asymptotically close to map fll.ip . Thus, 
using the results of [T] (see also [TT] for more generic statement), we deduce the birth of 
discrete Lorenz attractors in close systems. 

In the case of homoclinic tangencies to the saddle hxed point O of a three-dimensional 
diheomorphism T the birth of Lorenz attractors was proved in the cases when; 

1) |43| . the point O is a saddle-focus with the unit Jacobian (saddle-focus of conser¬ 
vative type). 

2) H, the hxed point is a saddle with the unit Jacobian and the quadratic tangency 
is non-simpl^. 

We note that the condition on Jacobians in all these cases is necessary for the existence 
of a non-trivial (three-dimensional) dynamics in the neighborhood of the homoclinic orbit 
|14) . Otherwise, if, for example, one has J < 1, all three-dimensional volumes will be 
contracted under the iterations of map T (near point O) and, hence, the dynamics of hrst 
return maps for large k will be ehectively two-dimensional, or even one-dimensional. 
Recall that, by dehnition [T3], the effective dimension de of a bifurcation problem equals 
n, if periodic orbits with n multipliers equal ±1 can appear at bifurcations but no orbits 
exist with more than n unit multipliers. 

Formally, if to consider three-dimensional diheomorphisms with homoclinic tangencies 

i.e. \\DTl^as II < La^ and | det DT^eu \ > Lv^ for some constants L > 0,0 < a < 1 < v and all positive k. 

^The definition of simple homoclinic tangency can be found in |12) . In particular, it assumes that the 
so-called extended unstable invariant manifold intersects transversely the leaf of the strong stable foliation 
in the point of tangency. The main cases of non-simple homoclinic tangencies in three-dimensional 
diffeomorphisms were considered in [2], see also condition D in §2 of the present paper. 
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to a hyperbolic saddle fixed point with \J\ = 1, then df, can be eqnal to 3 only in the 
following cases: (i) the point is a saddle-focus; (ii) the point is the saddle (all multipliers 
are real) and the tangency is not simple, and (hi) the point is a resonant saddle, i.e. it has 
two multipliers equal in the absolute value. Otherwise, the effective dimension is less than 
three since the direction of strong contraction is present for all nearby systems. 

Cases (i) and (ii) were considered in [13] and |1] respectively. In this paper we consider 
the new case (iii) when the saddle is resonant. Note that if the resonance Ai = A 2 takes 
place, we may perturb the map in such a way that the resulting map will have a saddle- 
focus hxed point with | J| = 1 and, hence, we can apply results of |13] to prove the birth 
of discrete Lorenz attractors. It not the case for the resonance Ai = —A 2 which is of 
independent interest. 

We consider the case when a hxed point O has multipliers A, —A, 7 such that 0 < A < 1, 
I 7 I > 1 and |A^ 7 | = 1. This means that O is a resonant saddle point of conservative type. 
Obviously, the bifurcation codimension of this problem is at least three and, as we will 
show, de = 3 in this case. 

We show that in the three-parametric families /^, fi = (/ii,/r 2 ,/^a) unfolding generally 
this type of a homoclinic tangency, in the parameter space there exist domains —)■ 
{/r = 0 } as fc —)■ cxD such that for fi G the hrst return map possesses the discrete 

Lorenz attractor. Recall that the map is constructed by the iterations of map /^, i.e. 
Tfc = /^, but the domain of its dehnition is a small box cTq near some homoclinic point. 

The paper consists of two paragraphs. In §2 we formulate our main result - Theorem [1] 
and construct the hrst return map of some small neighborhood of the homoclinic orbit. 
In §3 we prove Theorem [T] 


2 Statement of the problem and formulation of main 
results. 

We study bifurcations of three-dimensional diffeomorphisms of a special type (codimen¬ 
sion two) quadratic homoclinic tangency to a saddle hxed point with the unit Jacobian. 
Namely, we assume that the initial diheomorphism /o G > 5, satishes the following 
conditions; 


A) 


/ has a saddle hxed point O with real multipliers Ai,A 2,7 such that 0 < |Ai^ 2 | < 
1 < I 7 I and 

Jo = IA1A27I = 1 . 


B) The stable W^{0) and unstable W“(0) invariant manifolds of O have a quadratic 
tangency at the points of some homoclinic orbit Tq. 
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C) The saddle O is resonant in the sense that Ai = —A 2 = A > 0. 

Condition A means that the point O is a saddle of conservative type and dim W^{0) = 
2 and dim 1T“(0) = 1 . Condition C is an additional degeneracy of the saddle fixed point. 
We will consider smooth parameter families of diffeomorphisms (general unfoldings of 
conditions A-C), such that /o belongs to it for e = 0. 

Let U = 17(0 U To) be a sufficiently small fixed neighbourhood of Tq that is a union 
of a neighbourhood Uq of O and a number of neighbourhoods of those points of Tq which 
he outside Uq. Denote by Tq the restriction of the diffeomorphism onto Uq. We call Tq 
a local map. By a linear transformation of coordinates in Uq, map Tq can be written as 

(xi, 72 ,|/) = (AiXi, A2a;2,7|/) + h.o.t. 

The origin O = (0,0,0) is a fixed point of Tq, the stable manifold W^{0) is tangent at 
O to the (xi, 0 : 2 )-plane and the unstable manifold W“((9) is tangent at O to the ?/-axis. 
The intersection points of Tq with Uq belong to the set fl and accumulate to O 
at both forward and backward iterations. Thus, infinitely many points of Tq lie on 
and hT/oc- Let M+ e and M~ E be two such points and let M+ = fQ°{M~) for 
some positive integer uq. Let !!■*■ C Uq and 11“ C Uq be small neighbourhoods of points 
M~^ and M~ respectively. The map Ti = : 11“ 11+ is called a global map. 

From [32l |23 I2H1122] it is known that there exists a C^-change of coordinates (which 

is C"’“^-smooth in the parameters) bringing Tq to the so-called main normal form: 

xi = Ai(£)xi 0(||x|P|?/|) , 

X2 = X2{e)x2 + 0{\\xf\y\), ( 2 . 1 ) 

y = 7(^)1/+ 0(11x1111/12) _ 

The main peculiarity of this form is that in coordinates fl 2 .ip the stable and unstable 
manifolds of the saddle fixed point are locally straightened, their equations are : 
{y = 0}, hF“ : {xi = 0 ,X 2 = 0}. The main normal form also allows to obtain a quite 
simple representation of the iterations of Tq. The latter can be formulated as the following 
lemma: 

Lemma 1 j 3 ] For any positive integer k and for any sufficiently small e the map Tq (e:) : 
(2^0)?/o) {‘^kiVk) can he written in the following cross-form 

Xki - Af(e:)xoi = X^fki{xo,yk,e) , 

Xk2 - X2{e)xo2 = X’'fk2{xo,yk,e) , ( 2 . 2 ) 

yo - 'yie)-’^yk = f-’^T]kixo,yk,£) , 

where X and 7 are some constants such that A = A -|- 5, 7 = A| 7 “^| — 6 for some small 
(5 > 0 and functions f,k andrjk are uniformly hounded along with all derivatives up to order 

(r-2). 
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Next we construct the most appropriate form for the global map Ti for all small e. 
Let the chosen homoclinic points have coordinates M~^ = M~^{xf,X2,0) G and 

M~ = M~{0,0,y~) e where (xf)"^ + {x^Y ^ 0 and ?/“ > 0. At e = 0 we have 
that TiM~ = M+ and and are tangent quadratically at the point M+. 

Thus, the global map Ti at e = 0 can be written as the Taylor expansion near the point 

(xi = 0 ,X 2 = 0 , 1 / = y-): 

xi-xf = aiixi + 012X2 + bi{y - y-) + 0(||x|p) + 0(||x|||// - y-\) + 0 {{y - y-f) 
X2-xt = a2iXi + 022X2 + 62(2/ “ y~) + 0(j|xjn + 0(j|a;j| \y - y~\) + 0 {{y - y-f) 
y = CiXi + C2X2 + d(//-//“)^ + 0(||x|n + 0(||a;|||y-//-|) + 0 (|//-//-| 3 ) 

(2.3) 

The equation of curve Ti(lT;“^) at £ = 0 looks as follows (we put Xi = X 2 = 0 in fl2.3p h 

Xi-xt = hliy - y-) + 0{{y - y-f) 

X 2 -xt = b 2 {y - y-) + 0{{y - y-f) (2.4) 

y = d{y-y f+ 0{\y-y-\^) 


This is a parametric equation (with parameter [y — y~)) of the curve Ti(lT;“^) in a neigh¬ 
bourhood of . The equation of is ?/ = 0. Since the initial homoclinic tangency is 
quadratic, it follows that d 7 ^ 0, 6 ^ -|- 63 7 ^ 0. Moreover, map Ti(0) is a diffeomorphism, 
therefore 


Ji = det 


Oil ®12 bi 

O21 O22 &2 

Cl C2 0 


7^0 


(2.5) 


and, hence, c\ + 0 . 

At small £ the global map Ti{e) can be written in the following form (the Taylor 
expansion near the point {xi,X 2 ,y) = ( 0 , 0 ,//“(e)) 


xi-x+(e) = oiixi 012 X 2 + fei(//-//“(e)) + ... 

X 2 -x^(e) = O 21 X 1 -h 022 X 2 + 62 ( 2 /- l/■(^)) + ■ • • ( 2 . 6 ) 

y = //+(e) 4- CiXi -4 C 2 X 2 + d(y - y~{e)f + ... 


where x^( 0 ) = x^,xj( 0 ) = x^,//“( 0 ) = y~,y~^{0) = 0 ; all coefficients On ,... ,d depend 
(smoothly) on e; and we shift y~ into y~{s) in order to nullify the linear in y terms from 
the right side of the third equation. 

We assume that the following general condition holds 


61 C 162 C 2 7 ^ 0. (2.7) 

Its meaning is as follows: if at least one of these four coefficients is zero, then in any 
neighborhood of /o there exist maps having a nontransversal homoclinic orbit close to Tq 
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with a non-simple quadratic tangency (see mum for details); it composes an additional 
degeneracy which we do not consider here. 

Conditions A, B, C together with (12.71) dehne a codimension 3 bifurcation surfaces 
of diffeomorphisms with a quadratic honioclinic tangency. Hence, as a general unfolding 
we should consider a three-parameter families where the parameters fii, fi 2 and control 
the degeneracies imposed due to conditions B, A and C, respectively. 

Naturally, the splitting distance of manifolds W^{0) and with respect to the 

point is considered as the hrst governing parameter /ii. It is seen from (12. 6 p that 

hi = y^i^) ■ (2.8) 


The second parameter should control the Jacobian J 
Therefore, we dehne 


^2 = 1 - IA1A27I 


A 1 A 27 of ffj_ at saddle O^. 


(2.9) 


As the third parameter /is we consider the value that controls the difference between 
I All and IA 2 I, namely: 


h3 


|Ai(e)| 

|A2(5)| 


( 2 . 10 ) 


Thus, the family constructed above can be considered as a general unfolding 

of the corresponding homoclinic tangency to a resonant saddle, satisfying conditions A, 

B and C. 


Now we are able to construct the hrst return maps Tk using formulae (12. 2 p and (12.Op . 
As a result we will obtain a formula for Tk in the initial (small) variables (xi, X 2 ,y) G Uo 
and parameters /ii, /i 2 and Next, we rescale the initial variables and parameters 

{xi,X2,y) (Ai,A2,y) , (/ii,/i2,h3) {Ml, M 2 , M^) , 

with asymptotically small (as k 00 ) factors, in such a way that in the rescaled variables 
and parameters map Tk is rewritten as some three-dimensional quadratic map which 
contains asymptotically small (as k —)■ 00 ) terms. Moreover, new coordinates {Xi, X 2 ,Y) 
and parameters {Mi, M 2 , M^) can take arbitrary hnite values at large k (i.e. covering all 
values in the limit fc —)■ cx)). 

Our main result is the following theorem. 

Theorem 1 Let f 111 , 112,113 be the family under consideration. Then, in the (/ii,/i 2 ,/xa)- 
parameter space, there exist infinitely many regions accumulating at the origin as 
k ^ 00 such that the map Tk in appropriate rescaled coordinates and parameters is asymp¬ 
totically -close to the following limit map 
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where 


( 2 . 11 ) 


Xi = F, X 2 =^i, F = Ml + M2Xl + 5X2-F^ 


Ml = —+ AiCio:^ + A 2 C 2 XJ + o(A^)] (2.12) 

and 

M 2 = (b,ci + b2C2^ + •••), B = Ji(AiA27)"(1 + • • •) (2.13) 

3 Proof of Theorem [B 

Using fl2.6p and fl2.2p one can write the map = TiTq for snfficiently large k and small 
e in the form 


xi-xl = aii{X'^xi + X’^^kiix,y,e)) + ai 2 i>^ 2^2 + ^’'^k 2 ix,y,e))+ 

+bi{y - y~) + 0 {\y - y-\^ + X'^WxlWy - y-\) + X^^Wxf) , 

X2-x^ = ai2{X'lxi + X’^^kiix,y,e)) + a22{X2X2 +X’^^k2ix,y,e))+ 

+b2{y - y~) + 0 {\y - y-\^ + A^||x|||?/ - y-\) + X‘^^\\xf ) , ( 3 . 1 ) 

T^y-T''hk{x,y,e) = /ii+ ci(Ajxi+ A^^fci(x, ?/,£))+ 

+C2(A^X2 + X^ik2{,x, y, e)) + d{y - y-f+ 

+ 0 {\y — y~\^ + A^||x|||?/ — y~\ + A^^||x|p) . 

We shift coordinates ximw = xi - x+ + (plie), X2new = X2 - xj + (pUe), ynew = 
y — y~ + ' 0 fc(^); where ^’ki'^k = 0 {X^), in snch a way that the right sides of fl 3 .ip do not 
contain constant terms for the first two equations and linear in ynew terms for the third 
equation. Then fl 3 .ip takes the form 

xi = aiiAjxi + ai2A2X2 + biy + 0 {y‘^ + A^||x|| \y\ + A^||x|p) , 

X2 = a2iAiXi + a22A2X2 + 621/+ 0(|/2 + A''||x||||/| + A^||x|p) , 2) 

y-{ih)~''Tik{x + x-^+ (t)k,y^y~+ 'il^k,£) = Mk+ 

d'-y^y'^ + Ai7^CiXi + A 27 ^C 2 X 2 + 7^0(|?/|^ + A^||x|| \y\) + A^||x|p) 


where 

Mk = 7 ^[/ii + A 1 C 1 X++ A 2 C 2 XJ + o(A^)] (3.3) 
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Consider the third equation of 03.21) . First of all, we transform its left side. Namely, we 
write y - ( 7 / 7 )"''hfc = y + ( 7 / 7 )"^[ 7 ° + + vliv, y, e)] where 7 ^( 0 , e) = 0 , 

ril{0,e) = 0 and rjl = 0(||F^||). Next, we transfer constant term into the 

right side and join it to we substitute the value of x due to the first two equations 
of 03.2p into function if]\{x,e) and transfer the obtained expression into the right side. 
After this, all coefficients (in the third equation) get additions of order and a 

new linear term in ?/, p^y = /'y\~^)y, appears. By the shift of coordinates of the 

form {x^y) 1 —)■ {x,y) + 0 ([ 7 / 7 ]“^), we vanish both this linear term and constant terms in 
the right sides of the first and second equations. As the result, the left side of the third 
equation can be written as follows: 

y + {lh)~^0{y) + {^/-i)~^0(\\xy\\) = y{l + qu) + /-i)~^0(\\xy\\) 

where qk = ^([y/y]”^). After this, we can write system 03.21) in the form 

xi = aii\\xi +ai 2 Ax 2 + hiy + Oiy"^) + \^0{\\x\\\y\)+ \^0{\\x\\^) , 

X2 = a2i\\xi + a22>^2^2 + h2y + 0{y‘^) + \^0{\\x\\\y\) + \'^0{\\x\\^), 

s(l + ft) + (7/7)-‘0(|Sp) + (iH)-'‘ 0 (\\xv\\) = A.4+ 

+d 7 ^(l + Sk)y‘^ + ci\\^^xi + 02X21^x2 + Pkl^O{\\x\\‘^)+ 

+X^-i^O{\\x\\\y\) + -i^O{y^) , 


where Sk = 0(A^ + ly/y]”^), Pk = 0(X^ + ly/yl”^) and new Mk satisfies 03.3p . 

We perform a linear change of x variables to make zero the linear in y term in the 
second equation: 

_ ^2 _ _ 

^2new ^2 ~j ^Inew 7 Vnew ?/• 

Ol 

Then 03.4p is rewritten in the form 


= hiy + \’^0{\\x\\) + Oiy"^) , 

X2 = A2iX\xi + A 22 X 2 X 2 + 0{y^) + A^O(||x|||j/|) + A^^O(||a;|p) , 


2/(1 + gfc) + ( 7 / 7 ) ''0{\y\‘^ + \\xy\\) = + Afy^CiZ/fcXi + A^y'^CsXs + d 7 ''(l + 5^)//^+ 

Pk'-r''0{\\xf) + X^'j’^0{\\x\\\y\) + j’^oly^) , 


(3.5) 
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where 


^21 — ( 021 ~ Oil I + 1 O 22 ~ i—O 12 


h 


hi 


b_ 2 ^ 
hi Af 


^22 — 022 — 1 — 012 , 

Ol 

— / 1 _L "^2 
l^k — 1 + 


(3.6) 


61Cl 

Now we will vary Ai and A 2 in snch a way that the value of z/^ is asymptotically small 
as k ^ 00 . This is always possible via small changes of parameter because biCi 7 ^ 0 
and 62 C 2 7 ^ 0 and Ai = —A 2 in the initial moment. Then it is clear that 


^21 — 


A 

C 2&1 


O^Uk) 7 ^ 0 , 


where Ji is given by formula ( 12 .5p . 
Rescale the coordinates as follows 


y = 


7 (1 + gfc) 


Y ,xi = 


617 ’^{1 + qk) 


Xi , X2 = 


61^217 ^Xl{l + qk) 


d(l + Sfc) d{l + Sk) d{l + Sfc) 

Then system (13.51) is rewritten in the new coordinates as follows 

Xi = F + 0{X^) , 

X 2 = Xi + 0(7-"A-") , 

Y = Ml + M2X1 + BX2 -Y^ + 0 ( 7 -'=A-'^) , 


Xo 


(3.7) 


where formulas (I2.12p and (I2.13P are valid for Mi, M 2 and B. 

It is obvious that system (13.71) is asymptotically close to (12.lip when k —)■ 00 . □. 
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